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We theoretically study the optical absorption property of twisted bilayer graphenes with various 
stacking geometries, and demonstrate that the spectroscopic characteristics serve as a fingerprint 
to identify the rotation angle between two layers. We find that the absorption spectrum almost 
continuously evolves in changing the rotation angle, regardless of the lattice commensurability. The 
spectrum is characterized by series of peaks associated with the van Hove singularity, and the 
peak energies systematically shift with the rotation angle. We calculate the optical absorption in 
two different frameworks; the tight-binding model and the effective continuum model based on the 
Dirac equation. For small rotation angles less than 10°, the effective model well reproduces the low- 
energy band structure and the optical conductivity of the tight-binding model, and also explains 
the optical selection rule analytically in terms of the symmetry of the effective Hamiltonian. 



I. INTRODUCTION 

The recent advances in fabrication of atomically-thin 
materials have realized a new kind of two-dimensional 
superlattice in which the lattice mismatch between neigh- 
Atomic structure 



(a) = 21 .8° [(m,n)=(1 ,2)] (b) e = 5.09° [(m,n)=(6,7)] 




Brillouin zone 

(c) = 21 .8° [(m,n)=(1 ,2)] (d) 6 = 5.09° [(m,n)=(6,7)] 




FIG. 1: Atomic structures of TBGs with (a) = 21.8° and (b) 
= 5.09°. Dashed (orange) and solid (green) lines represent 
the lattices of layers 1 and 2, respectively. Brillouin zone of 
TBGs with (c) = 21.8° and (d) = 5.09°. Dashed (orange) 
and solid (green) large hexagons indicate the first Brillouin 
zone of layer 1 and 2, respectively, and thick small- hexagon 
is the folded Brillouin zone of TBG. 



boring layers gives rise to an additional potential mod- 
ulation. One example of such systems is twisted bi- 
layer graphene (TBG), in which two graphene layers are 
stacked with an arbitrary orientation. ^^^^^ In TBG, the 
interlayer interaction between two misoriented layers sig- 
nificantly modifies the low-energy band structure, arising 
novel electronic features distinct from intrinsic graphene. 
In decreasing the rotation angle, the interference be- 
tween two lattice periods produces a Moire pattern with 
a long wavelength, where the characteristic features such 
as band gaps and van Hove singularity appear in the far- 
infrared region, and the band velocity of Dirac cone is sig- 
nificantly r educed ^^iiiii — Recently, the band properties 
of TBG has been probed by Raman spectroscopy,^^^*" 
optical reflection spectroscopy,^^ angle-resolved photoe- 
mission spectroscopy,^^ and by terahertz time-domain 
spectroscopy.^^ 

The purpose of this paper is to reveal the optical ab- 
sorption properties of TBGs with various stacking ge- 
ometries. The optical absorption measurement is widely 
adopted for graphene-based systems to investigate the 
electronic structures. ^^^^ Theoretically, the optical ab- 
sorption for light incident perpendicular to the layer is 
related to the dynamical conductivity, which was calcu- 
lated for monolayer graphene, ^^^^ A^-stacked graphene 
bilayer and multilayers,^^^^ and also for TBG with a 
specific rotation angle j^^i^^ In this paper, we calculate 
the dynamical conductivity of TBG in a wide range of 
rotation angles, and demonstrate that the spectroscopic 
characteristics serve as a fingerprint to identify the stack- 
ing geometry. We find that the spectrum is characterized 
by series of peaks associated with the van Hove singu- 
larity, of which transition energy continuously shifts as 
the rotation angle is changed. Here the dynamical con- 
ductivity is calculated in two different frameworks; the 
tight-binding model and the effective continuum model 
based on the Dirac equation. For the latter, we develop a 
general treatment to derive the effective model from arbi- 
trary tight-binding parametrization. For small rotation 
angles less than 10°, the effective model nicely reproduces 
the low-energy band structure and the dynamical conduc- 
tivity of the tight-binding model, and also explains the 
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FIG. 2: Atomic structures of TBGs with (a) = 13.2°, (b) 11.0°, and (c) 9.43°. Dashed (red) and soUd (blue) parallelograms 
correspond to the Moire unit cell and rigorous superlattice unit cell, respectively. 



optical selection rule analytically in terms of the symme- 
try of the effective Hamiltonian. 



II. THEORETICAL METHODS 
A. Atomic structure and Brillouin zone 

Graphene is a single layer of carbon atoms arranged in 
a honeycomb lattice structure, of which unit cell includes 
two inequivalent sublattice sites, A and B. The stacking 
geometry of bilayer graphene is characterized by the rel- 
ative rotation angle combined with the lateral transla- 
tion S between the layers. Here we define the structure of 
TBG by rotating the layer 1 and 2 of AA-stacked bilayer 
around a common 5-site by — 6>/2 and +6>/2, respectively, 
and then translating the layer 2 relatively to the layer 1 
by S. We define ai = a(l,0) and a2 = a(l/2,>/3/2) as 
the lattice vectors of the AA-stacked bilayer before the 
rotation, where a ~ 0.246 nm is the lattice constant. The 
lattice vectors of the layer / after the rotation are given 
by a-^^ = R{'=f0/2)ai with ^ for / = 1,2, respectively, 
where R{0) represents the rotation by 0. 

When S is fixed to 0, the rotation ^ = and 60° give 
AA and AB stacking, respectively. 60° — ^ is equivalent 
to — ^ followed by a relative translation of the layer 2 
from A site to B site4^i^ Also, and —0 are mirror 
images sharing equivalent band structures. Therefore, it 
is reasonable to characterize the geometry of TBG by the 
combination of 6> (0 < 6> < 30°) and S. 

The lattice structure of TBG is not periodic in general 
angles because the periods of two graphene layers are 
generally incommensurate with each layer. But in some 
special angles where two periods happen to match, the 



structure becomes rigorously periodic giving a finite unit 
cell. This takes place when coincides with the angle 
between vi = mai + na2 and V2 = nai -f ma2 with 
certain integers m and n, because then the lattice points 
vi on the layer 1 and V2 on the layer 2 of the non-rotated 
bilayer graphene merge after the rotation 0/2 and —0/2^ 
respectively. The lattice vectors of the superlattice unit 
cell are thus given by^^ 
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and L2 = i?(7r/3)Li. The rotation angle is related to 
(m, n) by 
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The lattice constant L = |Li| = IL2I is 
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Figures [TJa) and [TJb) show the atomic structures 
of TBG with two different rotation angles = 21.8° 
[(m, n) = (1,2)] and 5.09° [(m, n) = (6,7)], respectively. 
Unless otherwise noted, ^ = is used. Figures [TJc) 
and md) show the corresponding Brillouin zone in the 
extended scheme. In each figure, dashed (orange) and 
solid (green) large hexagons correspond to the first Bril- 
louin zones of layer 1 and 2, respectively, and thick small 
hexagon to the folded Brillouin zone of TBG. i^^^) and 
K'^^"^ denote the two inequivalent valleys of layer I. The 
four valleys K^^\ K'^^\ K^'^\ and K'^'^^ of the two layers 
are folded back to the two Dirac points, K and K' ^ in 
the folded Brillouin zone^i^ 
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When the rotation angle is smah, the mismatch be- 
tween the lattice vectors of the two layers gives rise to a 
Moire pattern with a long spatial period as seen in Fig. 
[TJb).^^'^^ The local lattice structure near a certain point 
r approximates a non-rotated bilayer graphene with dis- 
placement which depends on the position as 



S{r) = 2sin((9/2)(e^ x r), 



(4) 



where r is measured from the center of the rotation, and 
Bz is a unit vector perpendicular to the plane. The period 
of the Moire pattern can be obtained by the condition 
that ^(L-^) coincides with a primitive lattice vector of the 
original AA-stacked bilayer. We may choose 
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The Moire superlattice vectors can be always de- 
fined for any 0^ even when the lattice structure is in- 
commensurate. At commensurate angles, the rigorous 
superlattice period L is exactly \m — n\ times bigger than 
the Moire period Lm- In Fig. [21 we illustrate the lat- 
tice structures of l9 = 13.2° [(m, n) = (2,3)], = 11.0° 
[(m,n) = (5,7)] and 6 = 9.43° [(m, n) = (3,4)]. We 
can see that the atomic structure exactly matches the 
Moire pattern in 13.2° and 9.43° (|m — n\ = 1), while in 
11.0°, the exact period L is twice as large as Lm since 
\m — n\ = 2, and accordingly the atomic structure is 
slightly different between neighboring units in the Moire 
pattern. 



B. Tight-binding model 

We calculate the eigenenergies and eigenfunctions in 
TBG using the tight-binding model for pz atomic or- 
bitals. The Hamiltonian is written as 



= - ^ t(Ri - R^-)|Ri)(Rjl + H.c, 



(7) 



where and |Ri) represent the lattice point and the 
atomic state at site z, respectively, and t(R^ — Rj) is the 
transfer integral between the sites i and j. We adopt an 
approximation 
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where ao = a/v^ ~ 0.142nm is the distance of neighbor- 
ing A and B sites on monolayer, and ~ 0.335 nm is 
the interlayer spacing. V^^^ is the transfer integral be- 
tween the nearest-neighbor atoms of monolayer graphene 
and V^p^ is that between vertically located atoms on 
the neighboring layers. Here we take V^^^ ^ — 2.7eV, 
Vpp^ ~ 0.48 eV, to fit the dispersions of monolayer 
graphene and A^-stacked bilayer graphene^ii^ is the 
decay length of the transfer integral, and is chosen as 
0.184a so that the next nearest intralayer coupling be- 



comes 0.11/; 



14,50 



The transfer integral for d > 4ao is 



exponentially small and can be safely neglected. 

C. Effective continuum model 

When the rotation angle is small and the Moire super- 
lattice period is much larger than the lattice constant, the 
interaction between the two graphene layers has only the 
long- wavelength components, allowing one to treat the 
problem in the effective continuum model. The contin- 
uum approaches for TBG have been introduced in several 
literatures Jiii^ii^i^ Here we develop a general treatment 
to construct an effective model directly from the tight- 
binding Hamiltonian in Eq. ([7]). To construct the Hamil- 
tonian matrix, we define the Bloch wave basis of a single 
layer as 

|k,B,) = -^^e*-^-<|RB,), (9) 

where the position Ra^(RbJ runs over all A{B) sites 
on the layer /(= 1,2), TV is the number of monolayer's 
unit cell in the whole system, and k is two-dimensional 
Bloch wave vector defined in the first Brillouin zone of 
monolayer on the layer I. 

Intralayer matrix element of each layer occurs only 
within the same wave vector, and it is given by 

/iA,A,(k) = {k,Ai\H\k,Ai) = h{k,0), 
/iA,B,(k) = {k,Ai\H\k,Bi) = hik,Ti), 

HbM^) = hA.AA^), (10) 

where ri = (2a2 — ai)/3 is a vector connecting from B 
site to A site, and 

/i(k, t)= ^ -t(niai + n2a2 + r) 

ni ,722 

X exp [ik • (niai + n2a2 + r)] . (11) 

The low-energy spectrum of the monolayer graphene is 
approximated by effective Dirac cones centered at K and 
points.^^^'^ We take K = (27r/a)(-2/3, 0) and = 
(27r/a)(2/3, 0) as the i^-points of non-rotated graphene. 
The K points of the layer / are then given by K^^^ = 
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R{tO/2)K and K'^^) = R{tO/2)K', with ^ for / = 1 
and 2, respectively. When k is close to either of K or K' 
the intralayer matrix element is approximately written 



-hv[{k, - K^^) - i{ky - K'y'^)]e-^^'^' (k ^ K), 

-hv[-{k, - K^'^) - i(ky - K'j'^)]e'^^'^ (k ^ KO, reciprocal lattice vector satisfying Lf • Gf = 2^5, 



(0^ 



where X and X' are either of A or 5, are the 

interlayer coupling in non-rotational bilayer in Eq. ([H]), 
and VLm = |Li^ x L2^| is the Moire superlattice unit cell. 
The derivation of Eq. ([16]) is detailed in Appendix [Al 
Ux^Xi [^^^(i*)] is periodic in r with the Moire superlat- 
tice periods, and therefore the matrix element Eq. (fT6|) is 
nonzero only when k' — k = niGj^ 



n2Gf where Of is 



(12) 



where 77^^^ = ±^/2 for / = 1 and 2, respectively. In the 
following we neglect the phase factor e~*^^'^ assuming 
6> <C 1. The parameter v is the band velocity of the 
Dirac cone, which is given in the present tight-binding 
parametrization as 
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where the first and the second terms in the bracket orig- 
inate from the hopping between the first and the second 
nearest AB pairs, respectively. The diagonal matrix ele- 
ment hAiAi (k) = hsiBi (k) are shown to be of the order of 
|k — Kp and |k — K'p near K and points, respectively, 
and will be neglected in the following. 

For the interlayer coupling, we first consider a non- 
rotated bilayer graphene with ^ = and a fixed lattice 
displacement 6. The unit cell is spanned by monolayer's 
lattice vectors, ai and a2, which are now shared by both 
layers. As the system has the same periodicity as the 
monolayer, the interlayer coupling occurs within states 
belonging to the same k. The interlayer matrix element 
is written as 

Ua,aA^,S) = {KA2\H\k,Ai)=u{k,S), 
Ub,b, (k, .5) = (k, B2\H\k, Bi) = u(k, ,5), 
Ub,a, (k, .5) = (k, B2\H\k, Al) = u(k, <5 - n), 
Ua,b, (k, .5) = (k, A2\H\k, Bi) = «(k, .5 + ti),(14) 

where 

ix(k, ^) = ^ -t(niai + n2a2 + d^ez + 5) 

ni,n2 

X exp [-ik ■ (niai + n2a2 + S)] . (15) 

i^(k, ^) can be immediately calculated by taking a sum- 
mation over some small n^'s, since t{d) rapidly vanishes 
in |d| > a. 

When is slightly shifted from zero to a small finite an- 
gle, the local lattice structure is approximately viewed as 
a non-rotated bilayer graphene, where the displacement 
S slowly depends on the position r in accordance with 
Eq. (|4j). Then the interlayer interaction couples wave 
vectors k and k' which are close to each other such that 
|k' — k| <C 2ix ja. The interlayer matrix element is ap- 
proximately written as 



{k',X^|F|k,Xi) 



k + k' 



-i(k'-k)-r 



,(16) 



and Tii is an integer. 

In TBG, the low-energy physics is still dominated by 
the states near K and K' points because the interlayer 
coupling is much smaller than intralayer coupling. Be- 
sides, the states near K and those near K' are far apart 
in the wave space when d is small, so that they not hy- 
bridized by the interlayer coupling. Therefore, we may 
consider two valleys separately in constructing the Hamil- 
tonian, and the factor (k + k')/2 in Eq. (p!6|) can be re- 
placed with K or K^ In the real space representation, 
the effective Hamiltonian near K is concisely written in 
the basis of IA2), 1^2)} as 
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V H2 
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with 



Hi 
U = 



u{K,S) u{K,S-Ti) 
u{K,6 + Ti) u{K,S) 



(18) 



where k = —id/dr^ a = {(Jx^cFy) is the Pauli matrices, 
AK^^) = K(^) - K, and 5 = ^(r) is defined in Eq. (j4]). 
When deriving Hi in Eq. ([18]) from Eq. ([T2]) , we replace k 
with k+K, i.e., measure the wavenumber relatively to the 
common point K for both layers. The Hamiltonian for 
K' is obtained by replacing K with and kx with —kx 
above. In the present choice of the tight-binding param- 
eters Vpp^^ Vpp^ and (5o, the effective interlayer coupling 
ii(K,^(r)) is approximately written in terms of only a 
few Fourier components as 



u{K,S{t)) ^ (O.lOSeV) X 1 + e 



Gi\ 
9 



(19) 



ii(K', S) is given by iz(K, ^)*. The expression of i^(K, 5) 
explicitly depends on the choice of K vector out of three 
equivalent corners in the Brillouin zone. 

In the /c-space representation, the Hamiltonian matrix 
can be written in the space of the single-layer bases at 
discrete /c-points k = ko + niGY + n2G2^, where ko is a 
vector defined in the superlattice Brillouin zone spanned 
by and G^^. ko = corresponds to the M point. 
To obtain the energy spectrum and eigen wave func- 
tion, we choose /c-points satisfying hv\]<.\<Emax with a 
sufficiently large £^max, and diagonalize the Hamiltonian 
within the limited wave space. To avoid a discrete change 
in the number of bases in varying ko , we adopt a soft cut- 
off which gradually reduces the matrix elements for the 
single-layer bases beyond ^max- 
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It is straightforward to show the Hamiltonian Eq. ((TTj) 
has a certain symmetry expressed as 

S-ij?effS = -H:^,, (20) 

where 

^=(-°/o-)- pi) 

This immediately concludes that if ip is an eigenstate 
of i^eff belonging to energy S?/^* is an eigenstate of 
energy —E. 

D. Dynamical conductivity 




X'(2)X'^1> K^2) j^(l) 



Using the eigen wave functions obtained by the tight- 
binding model or the effective continuum model, we cal- 
culate the dynamical conductivity 
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where the sum is over all states, S is the area of the sys- 
tem, f{e) is the Fermi distribution function, Sa {sfd) and 
I a) (1/3)) represent the eigenenergy and the eigenstate of 
the system, Vx = —{i/h)[x,H] is the velocity operator, 
and T] is the phenomenological broadening which is set to 
3meV in the following calculations. The optical absorp- 
tion intensity at photon energies hw is related to the real 
part of a{u). The transmission of light incident perpen- 
dicular to two-dimensional system is given by-^ 



27r , , 
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III. RESULTS AND DISCUSSIONS 



A. Band structure 



FIG. 3: Schematic band structures of TBG in the extended 
zone scheme along the line K'^^^-K'^^^-K^^^-K^^\ Symbols 
(i), (ii) and (iii) indicate the saddle points. 
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The band dispersion of TBG can be intuitively under- 
stood in terms of coupled four Dirac cones centered at 
the valleys K^^^ and K'^^^ {I = 1,2). Figure [3] shows the 
schematic band structures of TBG in the extended Bril- 
louin zone along the line K^^'^^-K^'^'^^-K'^'^^-K'^'^\ where 
Dashed and solid dispersion represent the energy bands 
of the layer 1 and 2, respectively. The interlayer coupling 
gives rise to band anticrossing at the intersection, and the 
resultant energy band is characterized by saddle points 
accompanied by the van Hove singularity in the density 
of states. The saddle points are classified into (i), (ii) 
and (iii) as shown in Fig. [3l where (i) and (ii) result from 
the band intersection of different layers, while (iii) orig- 
inates from the original monolayer's band structure. In 
the limit of (9 ^ 0, i^^^^ and i^^^^ (equivalently, K'^^^ and 
j^^(2) j closer, so that (i) approaches the Dirac points, 
while (ii) goes to the midpoint between K and K' in the 
high energy region, (iii) remains at the constant energy. 



FIG. 4: (a) Band structure and (b) density of states of TBG 
with = 3.89°. Solid (blue) and dashed (pink) arrows repre- 
sent the excitation corresponding to the major peaks in the 
optical absorption, and green dotted arrow is a process opti- 
cally forbidden (see text). 



Figures HJ^a) and HJ^b) show the band structure and 
the density of state (DOS), respectively, of TBG with 
6 = 3.89° [(m,n) = (8,9)] actually calculated by the 
tight-binding model Eq. ([7]). The band structure can be 
viewed as monolayer's Dirac cone folded into the super- 
lattice Brillouin zone with some band anticrossing at the 
zone corner i^Jiiiii^i^ — The lowest band exhibits a linear 
dispersion near K and K' ^ while we see a large splitting 
near the M point in energy from ±0.2eV to ±0.4eV, re- 
spectively, corresponding to the band anticrossing (i) in 
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FIG. 5: Atomic structures (top) and band structures (bottom) of TBGs with (a) = 9.43°, (b) 3.89°, (c) 2.65°, and (d) 1.47°, 
calculated by the tight-binding model (solid black lines) and the effective continuum model (dashed red lines). Dirac point 
energy is set to zero. 



Fig. The DOS has a sharp peak around ±0.2 eV, 

which is associated with the saddle point near the M 
point. 

Each energy band can be classified into either of those 
originating from monolayer's K region (i.e., K^^^ and 
or those from K' region (i.e., K'^^^ and K'^'^^), 
because two valleys are hardly mixed by the interlayer 
interaction in this small rotation angle. Here K and K' 
should not be confused with K and K' for the folded 
Brillouin zone. The monolayer's band near K and that 
near K' are independently folded into the same Brillouin 
zone without mixing with each other. Indeed, the low- 
est band in Fig. HJ^a) is composed of nearly degenerate 
branches, where dashed (pink) and solid (blue) lines are 
the bands from K and K' ^ respectively. These two bands 
are degenerate along K — T and M — K\ reflecting the 
C2 symmetry in the real-space lattice structure. 

Figure [5] shows the energy bands of TBGs with differ- 
ent rotation angles from = 9.43° down to 1.47°. The 
structures are similar to each other while the overall en- 
ergy scale shrinks as the rotation angle decreases, roughly 
in proportion to the size of the folded Brillouin zone. The 
width of splitting at the M point is about 0.2eV in every 
case, which is of the order of the interlayer coupling Vpp^. 
In small rotation angles less than 2°, the energy scale of 
the folded Dirac cone becomes comparable to the band 
splitting so that the band velocity near the Dirac cone 



is significantly reduced from the monolayer ' s In 
Fig. [5l we also plot the band energies calculated by the 
effective continuum model Eq. (fT7|) as dashed (red) lines, 
to be compared with solid (black) curves obtained by 
the original tight-binding model. We see the low-energy 
band structure agrees quite well, except that the effective 
model fails to reproduce a small electron hole asymme- 
try in the original model, since it assumes the symmetric 
Dirac cone for the intralayer Hamiltonian. 



Figure [6] shows the DOS of TBGs with various rotation 
angles 0° < ^ < 30° with (a) wide and (b) narrow energy 
ranges. To each curve, we append the DOS of uncoupled 
bilayer graphene (i.e., twice of monolayer's) as a light 
gray line. We observe a number of characteristic peaks 
associated with van Hove singularities (i), (ii) and (iii) 
argued in Fig. [3l As the rotation angle increases, the 
peaks (i) move away from the Dirac points, and the peaks 
(ii) move towards the Dirac points, while the peaks (iii) 
stay at almost constant energy. In small rotation angles, 
we see a number of additional peaks since the interlayer 
coupling of the higher order in the Moire wave number 
becomes significant. 
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FIG. 6: DOS of TBGs with various rotation angles 0° < ^ < 30° in (a) wide and (b) narrow ranges of the energy. The DOS 
of uncoupled bilayer graphene (i.e., twice of the monolayer's DOS) is shown as light gray lines. Peaks marked with symbols 
correspond to three different types of van Hove singularity (see text). 



B. Optical absorption 

In Fig. [3 we plot the dynamical conductivity of TBG 
with 6 = 3.89° calculated by the tight-binding model. 
The conductivity is plotted in units of 



16 h 



(24) 



which is the universal dynamical conductivity of mono- 
layer graphene at linear band regime, where gs = and 
Qy = 2 are the spin and valley {K^ K') degeneracy, 
respectively. ^^'^^^^ The spectrum is characterized by a 
peak around 0.5 eV, and considerable reduction right be- 
low the peak energy. Otherwise the conductivity is close 
to 2(Jmono- The top and bottom panels numbered from 
(1) to (4) show spectral weight maps at specific photon 
energies, which highlight the wave vectors that contribute 
to the optical transition. A sudden rise in the conduc- 
tivity between (3) and (4) is due to bright spots near M 
point in the weight map. This actually corresponds to 



the transition from the saddle point in the lowest valence 
band to the second conduction band, which is marked 
with a solid (blue) arrow in the band diagram in Fig.lU^a). 
The similar transition indicated by dashed (pink) arrow 
occurs at slightly higher energy due to the electron-hole 
asymmetry. Below the peak energy, the dynamical con- 
ductivity is significantly reduced and becomes lower than 
2crmono7 because the number of available states largely 
decreases due to the band anticrossing. 

It should be noted that the transition does not occur 
between the saddle points of the lowest conduction and 
valence bands [dotted (green) arrow in Fig.HJ^a)], because 
it is optically forbidden in this particular system. This 
can be clearly explained by the effective continuum model 
as following. The lowest electron and hole states at M 
point, which are connected by dotted arrow in Fig. 21 
are an electron-hole pair related by the symmetry of Eq. 
(|2Q|) , and thus these wavefunctions are written as and 
Tjip* . The matrix element of Vx between the two states 
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terized by a single peak belonging to group (i) similarly 
to Fig. [71 and its transition energy monotonically shifts 
with the rotation angle. In ^ = 1.47°, the spectrum 
exhibits a complicated structure in accordance with the 
strong band deformation observed in Fig. O In Fig. [D^b), 
we also present the spectrum of the effective continuum 
model as dashed (red) curves, to be compared with the 
original tight-binding calculation. The results agree quite 
well except for the peak splitting due to the electron-hole 
asymmetry, which is pronounced in ^>10°. 

As argued in Sec. Ill Al the rigorous superlattice period 
L in Eq. (|3]) discontinuously changes depending on the 
commensur ability of the lattice periods. Nevertheless, 
the DOS in Fig.[6]and the optical absorption spectrum in 
Fig. [8] show that the peak structure almost continuously 
evolves with the rotation angle, suggesting that the exact 
lattice commensurability is not quite important for the 
physical property. Although we cannot rigorously handle 
incommensurate TBGs due to infinite unit cell size, its 
optical spectrum should be approximated by interpolat- 
ing those of commensurate TBGs with similar rotation 
angles. The only property in which the actual lattice pe- 
riod L matters is found as tiny peaks in the conductivity 
indicated as dashed (blue) circles in Fig. [51 which exist 
only when \m — n\ > 1. They are related to the transition 
at the corner of the exact superlattice Brillouin zone. 



FIG. 7: (Middle panel) Dynamical conductivity of TBG with 
= 3.89° as a function of the transition energy. Arrows 
indicate the excitation energies of the transitions shown in 
Fig. mja). (Top and bottom panels) Spectral weight maps in 
the superlattice Brillouin zone at several transition energies. 



is obviously zero, because 




Mz) = 0, 



(25) 



and thus the transition is optically inactive. This sym- 
metry does not limit the optical selection rule at other 
points than M [the origin of the wave number in the ef- 
fective Hamiltonian Eq. ([T7|) ]. since and St/;* generally 
reside at different Bloch wave vectors and are not con- 
nected by the optical transition. 

In Fig. [Hl^a), we plot the optical absorption spectra at 
various rotation angles in a wide frequency range. The 
spectrum exhibits characteristic conductivity peaks rang- 
ing from terahertz to ultraviolet frequencies. The peaks 
are again classified into three groups similarly to the den- 
sity of states. When the rotation angle increases from 0° 
to 30°, the peak (i) [(ii)] moves to higher (lower) en- 
ergies, while the peak (iii) remains unchanged. Figure 
[8jb) shows magnified plots of the low-frequency range 
for several small angles. There the spectrum is charac- 



C. Effect of lattice displacement 

The electronic structure depends on the lattice dis- 
placement particularly when the superlattice period 
L is comparable to the atomic scale. This effect is 
maximum at 6> = 0, where the lattice constant L co- 
incides with a, and the A A stacking is transformed to 
the AB stacking by a translation S. The TBG with 
(m^n) = (1,2) {0 = 21.8°) has the next smallest primi- 
tive unit cell of L = ^/7a. Here we focus on two distinct 



TBGs of (m,n) = (1,2) with ^ = and (a^ 



(2) 



)/3, 



shown in Figs. [9]^a) and (b), respectively, to investigate 
the effect of the relative translation on the optical ab- 
sorption. These two structures can be also described by 
two different rotation angles 6 and 60° — 6. Figures [9fc) 
and (d) show the electronic structure and the dynami- 
cal conductivity, respectively, of the two different TBGs. 
We observe that the energy bands are almost equivalent, 
while there is a tiny difference in the structure near M 
point of which the energy scale is a few tens of meV. In 
the dynamical conductivity, this is reflected in a small 
difference in the peak structure near the energy of 3 eV. 
The inset in Fig. [9]^c) magnifies the band structure in 
the vicinity the Dirac points. We observe a difference of 
the energy scale about a few meV, where one takes the 
form like AA-bilayer and the other like AB-bilayer^i^i^ 
while it gives no noticeable difference in the dynamical 
conductivity in Fig. [9]^d). The dependence on the lat- 
tice displacement becomes even smaller as the superlat- 
tice period L becomes larger, and is completely absent in 
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FIG. 8: Dynamical conductivities of TBGs with various rotation angles in (a) wide and (b) narrow frequency ranges, calculated 
by the tight-binding model (solid black lines) and the effective continuum model (dashed red lines, only for (b)). Peaks marked 
with symbols represent the excitations associated with the van Hove singularity in Fig. \6\ Dashed (blue) circles indicate the 
tiny peaks which appears only when the actual lattice period L is larger than the Moire period Lm (see text). 



any incommensurate angles. Thus, concerning the energy 
range of interest, we conclude that the optical spectrum 
of TBG does not depends much on the lattice displace- 
ment, except for 6> = 0. 



low-energy band structure and the dynamical conductiv- 
ity of the tight-binding model for 6 < 10°, and it also 
explains the optical selection rule analytically in terms 
of the symmetry of the effective Hamiltonian. 



IV. CONCLUSION 



We theoretically investigated the optical absorption 
properties of TBGs with various stacking geometries us- 
ing the tight-binding model and the effective continuum 
model. We showed that the spectrum is characterized by 
series of absorption peaks associated with the van Hove 
singularities in the band structure, and the peak energies 
systematically shift in changing the rotation angle. The 
optical spectrum almost continuously evolves in changing 
the rotation angle regardless of the rigorous commensura- 
bility between two layers, suggesting that the optical ab- 
sorption measurement provides a convenient way to iden- 
tify the rotation angle of TBG. We developed the effec- 
tive continuum model based on the tight-binding model 
used here, and demonstrated that it well reproduces the 
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FIG. 9: Atomic structures of TBGs of 6* = 21.8° with different 
translation vectors, (a) ^ = and (b) (a^^'* + 8i^^)/3. (c) 
Band structures and (d) dynamical conductivities of the two 
distinct TBGs. 



Appendix A: Derivation of effective interlayer 
coupling 



Here we derive the effective interlayer matrix element 
Eq. (p!6|) for TBG with small rotation angles. The matrix 
element of the Hamiltonian Eq. ([7]) between the single 



layer bases on the different layers is explicitly written as 
{k',X^|F|k,Xi) 

= ^ -*(^^^ - RxJ exp[-ik' ■ Rx^ + ik ■ RxJ 



: exp [-ik ■ (R^^ - RxJ] exp 



-iAk 



where 



k = 



Ak = k' - k. 



(Al) 



(A2) 



Since the Moire lattice constant Lm is much larger than 
a in small rotation angles, we only need to consider k and 
k' which are close to each other, or |Ak| <C 2ix ja. Then 
the last exponential term in Eq. ()Aip slowly varies in the 
lattice position, while other terms change in the atomic 
length scale. To separate out the long wave component, 
we introduce a smoothing function ^(r) which satisfies 
the following conditions: ^(r) varies in r in an inter- 
mediate length scale which is much larger than the 
lattice constant a, but much smaller than the Moire su- 
perlattice period Lm- ^(r) is a peak centered at r = 0, 
and rapidly decays in r>/^. The area is normalized as 



/ 



M, 



(A3) 



where is the Moire superlattice unit cell, and the inte- 
gral is taken over the whole system area unless otherwise 
stated. Almost equivalent ly, we have 



(A4) 



where X is either of A\^B\^ A2 or and ^Iq = l^ii x 
is the area of monolayer's unit cell. 

Using Eq. ()A4p , the matrix element Eq. (jAip is written 

as 



(k',x;|i/|k,A'i) 

X exp [-ik • (Rx^ - RxJ] exp 



-iAk- 



t{Kxi - 



(A5) 



Now the argument (Rx^ + Rxi)/2 in the last can be 
replaced with R^i, because the last exponential term 
varies slowly with the length scale of Lm, and also the 
hopping integral t{Ilx^ — Rxi ) occurs only in the atomic 
scale distance. This is further replaced with r, since the 
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smoothing factor ^(r — RxJ works as the delta function where we used Eq. ()A4p and / dr/{Nflo) = Jq^ dr/flu- 

for the slowly varying function with the scale Lm- By The last equation is the final result of Eq. (pi)) , 
including this, we obtain 

{k',X^|F|k,Xi) 



X 



[-t(Rx^ - RxJ] exp [-ik ■ (R^f. - RxJ] 

R T^/ 

^ / dre-^^^-[/x^x, [k,^(r)], (A6) 
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